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' Aims. We examine the uncertainties in current planetary models and we quantify their impact on the planet cooling 

CO ' histories and mass-radius relationships. These uncertainties include (i) the differences between the various equations 

of state used to characterize the heavy material thermodynamical properties, (ii) the distribution of heavy elements 
within planetary interiors, (iii) their chemical composition and (iv) their thermal contribution to the planet evolution. 
Our models, which include a gaseous H/He envelope, are compared with models of solid, gasless Earth-like planets in 
order to examine the impact of a gaseous envelope on the cooling and the resulting radius. 

Methods. We find that for a fraction of heavy material larger than 20% of the planet mass, the distribution of the 
^ , heavy elements in the planet's interior affects substantially the evolution and thus the radius at a given age. For planets 

$— ( ' with large core mass fractions 50%), such as the Neptune-mass transiting planet GJ436b, the contribution of the 

C/3 I gravitational and thermal energy from the core to the planet cooling history is not negligible, yielding a ~ 10% effect 

, on the radius after 1 Gyr. 

Results. We show that the present mass and radius determinations of the massive planet Hat-P-2b require at least 200 
of heavy material in the interior, at the edge of what is currently predicted by the core-accretion model for planet 
formation. As an alternative avenue for massive planet formation, we suggest that this planet, and similarly HD 17156b, 
may have formed from collisions between one or several other massive planets. This would explain these planet unusual 
high density and high eccentricity. We show that if planets as massive as ~ 25 Mj can form, as predicted by improved 
core-accretion models, deuterium is able to burn in the H/He layers above the core, even for core masses as large as ~ 
100 M^. Such a result highlights the confusion provided by a definition of a planet based on the deuterium-burning 
limit. We provide extensive grids of planetary evolution models from 10 to 10 Mjup, with various fractions of heavy 
^SJ ■ elements. These models provide a reference to analyse the transit discoveries expected from the CoRoT and Kepler 

' missions and to infer the internal composition of these objects. 

QQ ' Conclusions. 
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1. Introduction rial such as rock or iccQ is not a surprise, since this is 

' 1 r 1 1 1 8- well known property of our own Solar System planets. 

The number of newly discovered exoplanets transitmg their Moreover, the presence of an icy/rocky core and of overso- 

parent star keeps increasing continuously, revealing a re- metallicity in the envelope is an expected consequence 

markable diversity m mean densities for planet masses ^^^^ ^.^^^^ accepted planet formation scenario, the 

ranging from^ Neptune masses to several Jupiter masses, ^o-called core- accretion model (Alibert et al. 2005a, and 

A large traction oi these transits exhibit a mean density references therein) 

significantly larger than that of an object essentially com- ' ^."^^^ this' expectation, it becomes mandatory to take 

posed of gaseous H/He, like brown dwarfs or stars, mdi- -^^^ ^^^^^^^ ^ ^^^^^^^ enrichment in planetary models 

catmg a composition substantially enriched m heavy ele- ^^^^^^^ ^^^^ ^^^^ identification of current and 

ments. The first compelhng evidence for such a significant forthcoming observations of extra-solar planets. Many ef- 

enrichment was provided by the discovery of a Saturn mass ^^^^^ ^^^^^^^ modeling of massive terrestrial 

planet, HD149026b, with such a small radius that 2/3 of ^^^^^^^^ essentially composed of soHd material (Valencia et 

the p anet s mass must be composed of elements heavier 2OO6; Sotin et al. 2007; Seager et al. 2007) and jovian 

than He (Sato et a - 2005) . Another remarkable discovery ^^^^^^ ^.^^ ^^^^ envelope and a substantial metal enrich- 
IS the case of GJ 436b, a ~ 22 Neptune-hke planet 

with a radius comparable to that of Uranus or Neptune 1 Under usual planet formation conditions, the word "rocks" 
(Gillon et al. 2007a). Such a radius implies an inner struc- refers primarily to sihcates (Mg-, Si- and O-rich compounds) 
ture composed by more than 90% of heavy elements. That whereas the term "ice" includes collectively H2O, CH4 and NH3, 
exoplanets can be substantially enriched in heavy mate- water being the most important of these three components. As 
will be discussed in §3, the term "ice" may be inappropriate in 



Send offprint requests to: I. Baraffe some cases, as water could be under a liquid or gaseous form. 
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ment (Baraffe et al. 2006; Guillot et al. 2006; Burrows et 
al. 2007; Fortney et al. 2007). Because of remaining uncer- 
tainties in the input physics describing the planetary struc- 
tures, and of many unknown quantities such as the total 
amount of heavy elements, their chemical composition and 
their distribution within the planet's interior, these mod- 
els are based on a number of assumptions and thus retain 
some, so far unquantified, uncertainties. The main goal of 
the present paper is to analyse and quantify these uncer- 
tainties and to explore as precisely as possible the impact of 
the heavy material contribution on the planet structure and 
evolution. We will focus on planets with mass larger than 
10 M^, the expected limit for the gravitational capture of 
a gaseous H/He envelope and atmosphere (Mizuno 1980; 
Stevenson 1982; Rafikov 2006; Ahbert et al. 2006), which 
has a significant impact on the evolution. Below this limit 
mass, the objects are essentially solid bodies with no or a 
teneous gaseous envelope and their mass-radius relation- 
ship has been studied recently by several authors (Valencia 
et al. 2006, 2007; Sotin et al. 2007; Seager et al. 2007). This 
study is motivated by the level of accuracy on planetary 
mass and radius measurements which is now reached with 
ground-based (HARPS, VLT) and space-based instruments 
(HST, SPITZER). Observations are expected to reach an 
unprecedented level of precision in the near-future with 
COROT, KEPLER, and on a longer term with GAIA. The 
latter project will measure distances and thus stellar radii 
with high precision, removing one of the main sources of 
uncertainty in planetary radius measurements. This race 
for precision is motivated by the possibility to infer with 
the best possible accuracy the inner composition of an exo- 
planet, with the aim to better understand planet formation 
and to identify the presence of astrobiologically important 
material such as liquid water. In this context, it is crucial 
to quantify the uncertainties in the structure and evolu- 
tion planetary models used to analyse these observations. 
In section §2 and §3, we examine the main input physics 
and assumptions used in structure and evolutionary models 
available in the literature. In section §4, we analyse quan- 
titatively the impact of these assumptions. We focus on 
specific cases such as HD149026b and GJ436b in §5 and on 
super Jupiter planets in §6. Our various planetary models, 
covering a wide mass range and including different levels of 
heavy element enrichments, are presented in §7. Discussion 
and perspectives follow in §8. 

2. Uncertainties and assumptions in the modelling 
of extra-solar planets 

2.1. Distribution of heavy material within the planet 

According to a recent study devoted to the structure of our 
giant planets (Saumon & Guillot 2004), Jupiter should have 
a total amount of heavy elements ranging from 8 to 39 M0 , 
i.e a metal mass fraction Z ^ 2.5% to 12%, with a maxi- 
mum core mass of 11 and a maximum envelope metal 
mass fraction Z^nv = 12%. For Saturn, the same study sug- 
gests a total mass of heavy elements ranging from 13 to 28 
Mq, i.e Z ^ 13% to 29%, with a maximum core mass of 
22 M0 and a maximum Zcnv ^ 8%. These properties can 
be understood within the standard general framework of 
giant planet formation : as a core is growing in mass due to 
accretion of planetesimals, it reaches a critical mass around 
^ 6-10 M0 above which gas accretion begins (Mizuno 1980; 



Stevenson 1982; Rafikov 2006). During this gas accretion 
phase, planetesimals are still accreted and are either de- 
stroyed in the gas envelope or are falling onto the core, 
leading to further increase of the core mass. The fate of 
these accreted planetesimals, disrupted in the envelope or 
accreted onto the dense core, depends on their size, an un- 
known parameter in current planet formation models, and 
on the envelope mass. This general picture thus predicts 
that planets should have a dense core of heavy material 
(rock, water/ice) of, at least, a few Earth masses and can 
show different levels of heavy element enrichment in their 
envelope, depending on the accretion history (amount of 
gas accreted, size of planetesimals, etc.). Current mod- 
els based on the core-accretion scenario are able to match 
the core mass and the envelope metal enrichment derived 
for Jupiter and Saturn (Alibert et al. 2005b), but predict 
in some cases much larger heavy element enrichments in 
the envelope, depending on the planet's mass, with values 
as large as Zcnv > 50% for Neptune-size planets 10-20 
Me) (Baraffe et al. 2006). 

Despite this widely accepted picture, current extra-solar 
planet models often simply assume that all heavy elements 
are located in the central core and that the envelope is 
either metal-free, Zg^v = 0, (Burrows et al. 2007; Fortney 
et al. 2007), or has a solar metallicitjO, Zcnv — Zq (Guillot 
et al. 2006). This simplification is based on the assumption 
that whether the heavy elements are located in the core or 
in the envelope should not affect the planet's evolution. The 
validity of such an assumption, however, has never been 
examined, as will be done in the present paper. 

2.2. Thermal contribution to the EOS 

An other simplification found in planet modelling is the use 
of temperature-independent EOS, assuming that the heavy 
element material is at zero-temperature or at a uniform, 
low temperature (Seager et al. 2007; Fortney et al. 2007). 
This is certainly a good assumption when examining the 
structure of terrestrial-like planets, composed essentially of 
solid (rocky/icy) material (Valencia et al. 2006; Sotin et al. 
2007). This assumption, however, is not necessarily valid for 
the evolution of more massive planets. For these objects, the 
thermal and gravitational energy contributions of the core 
to the planet's cooling history are usually ignored and it is 
important to examine the impact of such a simplification. 

2.3. Heat transport 

Finally, a conventional assumption in planet modelling is 
to assume that the interiors (at least the gaseous envelope) 
of giant planets are homogeneously mixed, due to the dom- 
inant and supposedly efficient transport mechanism pro- 
vided by large-scale convection. Accordingly, the internal 
temperature gradient is given by the adiabatic gradient, 
since superadiabaticity is negligible. The validity of this as- 
sumption, however, has been debated for the interior of our 
own jovian planets for already some time (Stevenson 1985) 
and has been questioned again more recently in the context 
of transiting extra-solar planets (Chabrier & Baraffe 2007). 
It should be kept in mind, however, that, even though the 

^ Note also that different envelope helium mass fractions have 
been assumed, with Y = 0.25 (Burrows et al. 2007) or y=0.28 
(Fortney et al. 2007). 
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planetary internal heat flux can only be carried out by con- 
vection (Hubbard 1968), the assumption of large-scale, fully 
adiabatic convection in planetary interiors has never been 
proven to be correct and even slightly inefficient convec- 
tion can have a major impact on the planet's structure and 
evolution (Chabrier & Baraffe 2007). 

3. Equations of state for heavy elements 

3.1. The case of water 

As quickly mentioned in the introduction, the generic term 
"ice" may often be inappropriate to describe the thermo- 
dynamic state of water under planetary conditions. Indeed, 
depending on the temperatures prevailing in the parent pro- 
toplanetary nebula at the initial location of the planet em- 
bryo, water could possibly be initially under the form of 
solid ice, but it could also be under the form of liquid or 
vapor. If initially solid, it may also melt or vaporize under 
the conditions prevailing in the planet's interior, or may 
also dissociate under the form of ionic melts at high pres- 
sures and temperatures (Schwegler et al. 2001). In fact, 
the phase diagram of heavy elements under the pressure 
and temperature conditions characteristic of the consid- 
ered planet interiors is largely unknown, as only part of 
this diagram is presently accessible to high-pressure exper- 
iments or computer numerical simulations. Water, for in- 
stance, the dominant component after H and He, is known 
to exhibit a complex phase diagram with many stable or 
metastable (amorphous) phases and several triple and crit- 
ical points, because of the high flexibility of the hydrogen 
bonding. The melting line of water at high pressure and 
high temperature has been probed experimentally up to 
P = 3.5 X 10" dyn cm'^ (0.35 Mbar) and T = 1040 K. 
H2O has been found to remain solid at larger pressures and 
temperatures, suggesting that the melting line increases (in 
T{P)) at higher temperatures and extends up to at least 
T>1600 K (Lin et al. 2005). Depending on its extrapo- 
lation at higher temperatures, this line may intersect the 
internal profiles of some of our planet models, at least in 
the Super-Earth and Neptune mass range, so that both 
liquid and solid H2O might be present at some depth in 
the interiors of these objects, as predicted for ocean-planets 
(Selsis et al. 2007), while supercritical H2O is more likely 
to be liquid or gaseous in the hotter interiors of Saturn- 
like or larger planets. As mentioned above, water is also 
found from shock- wave experiments and first-principles cal- 
culations to dissociate into H-^O'^ -\- 0H~ ion pairs above 
~ 2000 K at 0.3 Mbar (Schwegler et al. 2001). In fact, the 
distinction between "gas", "ice" and "rock" may become 
meaningless under extreme conditions. In that case, what 
matters is the global amount of heavy elements. The term 
"ice", or "water", may thus more generically refer to the 
volatile forms of O, C and N, while "rock" refers primar- 
ily to silicates (Mg-, Si- and O-rich compounds) and "iron" 
refers to the rest (metal, oxide, sulfide or substituting for 
Mg in the silicates) (Stevenson 1985). 

3.2. EOS available in the literature 

For Earth-like planets, models can test a variety of com- 
plex heavy element compositions inspired by the knowl- 
edge of the structure of our own Earth (Valencia et al. 
2006, 2007; Sotin et al. 2007). Moreover, the EOS of ma- 



terials which may be found in planetary interiors (water, 
iron, dunite or olivine, etc..) are reasonably well deter- 
mined at zero-temperature or at 300 K (see Seager et al. 
2007). Unfortunately, for larger planets, the exploration of 
internal compositions is restricted to a few materials for 
which EOS are available and cover a large enough range 
of pressures and temperatures. The two most widely used 
EOS in this context (Saumon & Guillot 2004; Baraffe et 
al. 2006; Fortney et al. 2007; Burrows et al. 2007) are 
ANEOS (Thompson & Lauson 1972) and SESAME (Lyon 
& Johnson 1992), which describe the thermodynamical 
properties of water, rocks (olivine or dunite, i.e Mg2Si04, 
in ANEOS; a mixture of silicates and other heavy elements 
called "drysand" in SESAME) and iron. Figure [T] shows a 
comparison between these two EOS for the three aforemen- 
tioned materials for two temperatures, namely T ~ 300 K 
and T — 6000 K. Comparison is also shown with the zero- 
temperature EOS presented in Seager et al. (2007) for wa- 
ter, perovskite (MgSiOs) and iron. These authors use fits 
to experimental data at low pressure {P < 2 Mbar) and 
an improved Thomas-Fermi-Dirac theory at high pressure 
(P > 100 Mbar), where the contribution of the degenerate 
electron fluid becomes dominant. As stressed by these au- 
thors, the main difficulty is to bridge the gap in the pressure 
regime 2 Mbar^P^ 100 Mbar (2 x lO^^ _ ^q" dyne cm'^), 
which is the most relevant for planetary interiors. Their 
model EOS represents some improvement upon ANEOS 
and SESAME, in particular for water, but is valid only at 
zero-temperature. 

The presently used EOS models (ANEOS and 
SESAME), and thus the inferred planet internal structures, 
thus retain a significant degree of uncertainty in the ex- 
perimentally unexplored high-P and high-T domains. In 
these regimes, both EOS models, as well as the one used by 
Seager et al. (2007), are based on interpolations between 
experimental data at low or moderate density/temperature 
and well-known asymptotic limits, in general the Thomas- 
Fermi or more accurate density-functional type models, in 
the very high density, fully ionized limit. One can only hope 
that these interpolations do not depart too much from re- 
ality, as might be the case, for instance, if first-order phase 
transitions, which imply density and entropy discontinu- 
ities, occur in the regions of interest. For the static prop- 
erties, this assumption is probably reasonable at high tem- 
peratures, but might be more questionable near melting 
lines. When addressing the transport properties, like e.g. 
the thermal diffusivity or kinematic viscosity, the results in 
the interpolated regime are definitely more doubtful. 

As shown in Figure [1] at low (room) temperature, the 
various EOS agree reasonably well. For water, the agree- 
ment lies within less than 5% at P=0.1 Mbar and within 
less than 16% at P = 100 Mbar. For "rocks", keeping in 
mind that this term refers to different compositions in the 
various EOS, the agreement is comparable, although the 
difference can reach 27% at P=0.1 Mbar. For iron, the three 
EOS agree well at low pressure but differences as large as 
20% can occur at P=100 Mbar between SESAME and 
Seager et al. (2007). Such cool temperature conditions, how- 
ever, are more relevant to Earth-like planets than to the 
ones explored in the present study. 

Figure[T]also shows, for ANEOS and SESAME, the vari- 
ation of P(/o) with temperature, for conditions more suit- 
able to our planetary interior conditions. The thermal con- 
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Fig. 1. EOS for water, rock and iron for two temperatures: T=300K for ANEOS (red solid line) and SESAME (blue 
dash-dotted curve); T=6000K for ANEOS (red short-dashed curve) and SESAME (blue dotted curve). The long-dashed 
(magenta) curves correspond to the zero-temperature EOS of Seager et al. (2007; their Table 3). The inset in the lower 
panel (water) shows a zoom of the P{p) relation for log P < 0.5 



tributions predicted by these EOS are significant at P < 1 
Mbar, with a ~ 40-60% difference in P{p) between the 
T = 300 K and T = 6000 K isotherms, for iron and water, 
respectively. The differences keep increasing significantly 
for T > 6000 K. It is instructive to quantify the impact of 
these thermal contributions on the evolution of our plan- 
ets, and to determine whether neglecting the temperature 
dependence of the EOS is consequential or not. 

Whereas P(p) is the relevant quantity for the structure, 
the relevant one for the evolution is the entropy. Figure 
[5] portrays the P- and T-dependence of the entropy, for 
three isotherms and isobars, for water, for the ANEOS and 
SESAME EOS, under conditions relevant to the planets 
of present interest. The values for the H/He fluid (EOS 
of Saumon et al. 1995; hereafter SCVH EOS) are also dis- 
played for comparison. We see that, whereas the EOS agree 
reasonably well under Jovian-planet conditions, as expected 
as they reach the asymptotic, high-P, high-T regime accu- 
rately described by Thomas-Fermi-Dirac or more accurate 
density-functional theories, the difference can be substan- 
tial for conditions characteristic of the Neptune-mass do- 
main. Indeed, for this latter case, most of the interior lies 
in the interpolated regime where guidance from either ex- 
periments or numerical simulations is lacking. These dif- 
ferences between the EOS, of course, are amplified for the 



quantities involving the derivatives of the entropy, such as 
the adiabatic gradient or the specific heat. 

3.3. Treatment of metal enrichment in the core and in the 
envelope 

Our goal is to examine the impact on the planetary models 
and the inferred mass-radius relationships due to uncer- 
tainties in (i) the distribution of heavy elements within the 
planet interiors, (ii) their dominant chemical composition, 
(iii) the different EOS describing their thermodynamical 
properties and (iv) their thermal contribution. To achieve 
this goal, we have implemented the ANEOS and SESAME 
EOS for water, rock and iron (or what is so-denominated) 
in our planetary evolution code. Both EOS provide all ther- 
modynamic quantities relevant to the evolution of planets, 
including internal energy, entropy (in ANEOS) or free en- 
ergy (in SESAME), and all relevant derivatives. ANEOS 
also provides Rosseland and conductive opacities. When 
the SESAME EOS is used, conductive opacities are calcu- 
lated according to Potekhin (1999). The evolutionary cal- 
culations, including the presence of a dense core, proceed as 
described in Baraffe et al. (2006): the structure equations 
are integrated from the center to the surface; at the core 
boundary, the EOS is switched to the one characteristic of 
the gaseous envelope, and the change in chemical composi- 
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Fig. 2. Pressure (in dyne cm ^) and temperature (in K) 
dependence of the entropy (in erg/g/s), as obtained for wa- 
ter for the ANEOS EOS (solid hne) and the SESAME EOS 
(dash-hne). The entropy for a H/He fluid, with Y = 0.275 
obtained with the SCVH EOS is shown for comparison 
(dash-dotted line). Upper panel: for each EOS, three isobars 
are displayed (the corresponding pressure is given on the 
right hand side of the figure). Lower panel: three isotherms 
are shown for each EOS (corresponding temperatures given 
on the right hand side of the figure). The T and P domains 
characteristic of Neptune-like and Jupiter-like planet inte- 
riors are indicated. 



tion yields a density jump, but continuity in pressure and 
temperature is enforced. 

To account for the thermodynamic contributions of 
heavy elements in the H/He envelope, we have tested two 
different procedures. (1) Following the method described in 
Chabrier et al. (1992) and used in Baraffe et al. (2006), 
we first mimic the presence of metals with mass fraction 
Z by an equivalent helium fraction ycquiv = Y + Z in the 
H/He SCVH EOS, with Y the real helium mass fraction in 
the envelope. As mentioned in Chabrier et al. (1992), this 
approximation is reasonable as long as the mass fractions 
Y and ycquiv are small compared with unity and both pz 
and pHe are large compared with pn, where pi denotes the 
mass density of the i-component at pressure P. (2) The sec- 
ond, more general approach to describe the EOS of a mix- 
ture of various species in the absence of a reliable theory is 
to apply exactly the additive volume law (hereafter AVL), 
which is exact in the ideal gas limit, without restriction on 
the species mass fractions and densities (see Fontaine et al. 
1977 and Saumon et al. 1995 for extensive discussions of 
the validity of the AVL). In this approach, the interactions 
between the three different fluids, namely hydrogen, helium 
and the heavy element component, are not taken into ac- 
count (but interactions between particles in each of these 



fluids are treated properly) and the EOS of the mixture 
is simply the mass-weighted interpolation of each species 
contribution at constant intensive variables, P and T, plus 
the ideal entropy of mixing for the entropy term. Within 
the ideal volume law, the mass density of the mixture of a 
H/He fluid with a helium mass fraction Y plus some heavy 
element material with mass fraction Z at pressure P and 
temperature T thus reads: 



1 



(1-Z) 



p{P,T) pH/He^p^T) p^{P,T) 



(1) 



The extensive variables, e.g. the internal energy and specific 
entropy, thus read: 

U{P,T) = (1 - Z)U"'"^{P,T) + ZU^{P,T) (2) 

S{P, T)^{l-Z) S^/^^iP, T) + Z S^{P, T) + 5„ii,(P, T) 

(3) 

where the EOS of the H/He component is given by the 
SCVH EOS while the one of the Z-component is described 
by either the ANEOS or the SESAME EOS. All along the 
present calculations, we have taken the cosmic helium frac- 
tion y=0.275. The details of the calculation of the con- 
tribution due to the ideal entropy of mixing are given in 
Appendix A. This term is found to contribute non negligi- 
bly to the total entropy S. Depending on the mixture and 
the P-T range, it can amount to 10%-20% of S. This is 
consequential when calculating the internal structure of a 
planet, whose interior is essentially isentropic, at a given 
age. 

For the evolution, the relevant quantities are the vari- 
ation of the entropy with time and its derivatives w.r.t P 
and r, which give the adiabatic gradient: 



_ d\ogS_ ,,d\ogS_ 
"""'^dlogP'^'^dlogT'' 



(4) 



Given the fact that the internal composition of the 
planet {NhjNhc, Nz, where Ni denotes the number of par- 
ticles of species i), remains constant along the evolution, 
the only P- and T-dependence of the mixing entropy term 
arises from the variation of the degree of ionization and 
thus of the abundance of free electrons. The degree of ion- 
isation of the heavy element component, however, is un- 
known, and its variation with temperature has been ignored 
in the present calculations. Therefore, only the variation 
with temperature of the number of electrons provided by 
H and He contribute to the variation with time of the mix- 
ing entropy. Within this limitation, the variation of Smix 
with P and T is found to be negligible for planets in the 
mass range of interest, for all levels of metal enrichment in 
the envelope, so that the entropy of mixing term does not 
contribute significantly to the evolution. 

4. Effect of the different treatments and 

distributions of heavy elements on the planet 
cooling history 

In this section, we analyse the impact of the localisation 
of the heavy elements within the planet on its structure 
and evolution. The heavy elements are distributed either 
in the core or in the gaseous envelope, and we consider 
in the present section mass fractions with Z < 50%. This 
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encompasses the level of enrichment of our giant planets 
Jupiter and Saturn and of previous theoretical studies de- 
voted to the analysis of transit planets (Guillot et al. 2006; 
Burrows et al. 2007). The effect of larger fractions of heavy 
elements {Z > 50%) will be explored in the next section, 
with the analysis of specific cases of planetary transits. To 
cover the largest possible range of conditions, we explore 
these effects on a Neptune mass (20 M^) and a Jovian 
mass planet (1 Mj = 318 M^). As a test case, we restrict 
the analysis in this section to planets irradiated by a Sun at 
0.045 AU, since the observational determination of the ra- 
dius is presently accessible only to transiting, short-period 
planets. We use outer boundary conditions derived from a 
grid of irradiated atmosphere models with solar metallic- 
ity (Barman et al. 2001). These models take into account 
the incident stellar flux, defined by Fine — ^i/^)^^*^ with 
Ri, and F* the radius and the flux of the parent star, re- 
spectively, and / the redistribution factor. The Barman 
et al. (2001) models use /=2, corresponding to the inci- 
dent flux being redistributed only over the dayside of the 
planet. Recent observations of the day- night contrast of ex- 
oplanets with Spitzer (see e.g. Knutson et al. 2007) and at- 
mopsheric circulation models (see the discussion in Marley 
et al. 2007), however, seem to favor a redistribution over 
the entire planet's surface, i.e f=l. For the specific case 
of HD209458b, Baraffe et al. (2003) found out that such a 
variation of Fine by a factor 2 has a significant effect on the 
outer atmospheric profile, but a small effect on the planet 
evolution. In a forthcoming paper, we will explore in more 
details the effects of the redistribution factor /, of different 
levels of irradiation and of heavy element enrichment in the 
atmosphere. 

4.1. Neptune-mass planets 

Figure [3] shows the evolution of the radius with time for 
a 20 planet with Z=50%, for different heavy element 
materials and for different distributions of these latter in 
the planet interior. We first consider models with a 10 M0 
core and Zcnv=0, i.e. all heavy elements are located in the 
central core of the planet. The results for pure water (solid 
line), rock (dashed line) and iron (dash-dotted line) cores 
are displayed in Fig. based on the ANEOS EOS for 
these materials. A comparison of these results with those 
obtained with the SESAME EOS for the same material 
(water and drysand) shows less than a 1% difference on 
R at a given time. For iron, the SESAME EOS does not 
provide the free energy and thus the entropy, the adiabatic 
gradient and other quantities required for the evolution. 
However, according to the comparison between the ANEOS 
and SESAME EOS portrayed in Fig. [TJ we do not expect 
significant differences on the evolution between the different 
iron EOS. 

For the aforementioned global heavy element enrich- 
ment, Z = 50%, the nature of the core material affects 
the cooling, and thus the radius evolution at the ^--^10% 
level. We will see in §5, however, that for larger heavy el- 
ement mass fractions, the impact of the core composition 
can be more severe. At 1 Gyr, the radius of the planet with 
a pure rocky (iron) core is smaller by 6% (11%) compared 
to the pure water core case (see Table 1) . For water and 
iron, the ANEOS EOS indicates whether the material is in 
a solid, liquid or melt phase. For dunite, this information 



is not provided. In all the models, water is always found to 
be in a liquid state. However, as mentioned in §3.1, most 
of our planet interiors lie in the extrapolated, supercriti- 
cal region of the phase diagram; the precise state of water 
in this regime depends on the extrapolation of the melt- 
ing curve and thus is presently undetermined. Moreover, as 
mentioned earlier, at high pressure and temperature, the 
meaning of "water" becomes loose and refers more gen- 
erally to protons and oxygen nuclei. The planet with an 
iron core undergoes a phase change from liquid to solid, 
according to ANEOS, at an age of ^ 1 Gyr. This transition 
starts in the most central part of the core, at temperatures 
T ~ 1.7 X 10"* K and p ~ 30 g.cm-^. 

Figure [SJa illustrates the effect of the localisation of the 
heavy elements within the planet's interior, with two lim- 
iting assumptions: (A) all heavy elements are in the core 
(Zonv=0) and (B) they are distributed all over the planet's 
interior (Mcore=0). In the latter case, we compare the two 
procedures mentioned in §3.2 to describe the thermody- 
namics properties of a mixture of H/He and heavy elements. 
Comparison between these two limiting cases should pro- 
vide the maximum effect on the planet's radius due to the 
unknown distribution of heavy elements within its interior, 
as long as large-scale adiabatic convection is considered as 
the heat transport mechanism. All the results portrayed in 
Fig. are done with water. We flnd the following results: 

— Under assumption (B), models based on a Yequiv, for 
such a high metal fraction, Z = 50%, lead to cooling 
sequences that differ drastically from the ones based on 
the AVL, with the SESAME EOS. This shows that the 
Yequiv simplification can not be used for such values of 
Z , as anticipated from the limitations of this assumption 
(see §3.3). 

— Case (A) and case (B), with the AVL and the SESAME 
EOS, yield similar coohng sequences, with ~ 4% differ- 
ence at t=l Gyr (7% at 5 Gyr). 

— In case B, we find a significant difference when using 
the AVL with ANEOS compared with all other cases 
30% difference in R for t > 1 Gyr compared with case 
(A), see Table 1). As mentioned in §3.2 and illustrated in 
Fig-El this stems from the different entropy dependence 
on (P, T) predicted by the different EOS in the relevant 
pressure regime indicated in the middle panel of Fig. [3) 
Although all the evolution sequences are forced to start 
from the same initial entropy state, the sequence based 
on the ANEOS EOS shows a much stronger variation 
of the entropy with time. This stems from the stronger 
variation of S with P and T predicted by this EOS, for 
the mass range characteristic of Neptune-like planets 
(see Fig. [2]). This is highlighted in the upper panel of Fig. 
m the sequence based on the AVL with ANEOS shows 
^ 10 times larger local gravothermal energy, —TdS/dt, 
at the beginning of the evolution. Consequently, this se- 
quence contracts and cools much faster than the other 
ones, loosing its internal entropy at a faster rate (lower 
panel of Fig. |4|) and reaching a significantly smaller ra- 
dius after 1 Gyr (see Fig.jSb). 

The reason why the difference between the entropy, 
S{P, T), obtained with the ANEOS and SESAME EOS 
(Fig. [2] and Fig. 3) is inconsequential when all the 
Z-material is located in the core, whereas it yields 
significantly different cooling sequences when the Z- 
component is mixed with the H/He envelope stems from 
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Table 1. Radius and central thermodynamic properties of a planet of mass Mp — 20 at 1 Gyr with a total 
heavy element mass fraction Z = {Mz /Mp)=50%. Results for different heavy element compositions (water=W, rock=R, 
iron=I) and EOS (aneos=a, sesame=s, SCVH=SC) and different heavy element distributions are given. 
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the fact that the evolution of a planet with a core is dom- 
inated by the entropy variation of the H/He envelope, 
\dS/dt\HHe ^ \dS/dt\z, a consequence of the smaller 
compressibility, [p{dP/dp)]~^ , and smaller specific heat 
of the Z-material compared with the H/He one because 
of its much larger mean molecular weight (see §4.3). 

Figure [S}; also shows the evolution of a model with a 
more realistic distribution of the heavy material in the in- 
terior, namely a 8 Mq core surrounded by a gaseous en- 
velope with a Zenv=16.6% heavy element mass fraction, 
which corresponds to a total heavy element enrichment of 
50%. First of all, we note that, for this value of Z, the cool- 
ing sequence base on the Ye/ / formalism (dash-dot line) is 
very similar (within < 2%) to the one based on the more 
rigorous AVL (long-dash line). This shows that the Yg// 
approximation for the treatment of a multispccies H/He/Z 
EOS can be used relatively safely up to Z « 20%. Second 
of all, sequences with Mcore = 8 M^, Zcnv=16.6% (with 
the AVL and the SESAME EOS or with a Yoquiv) and the 
ones with Mcore = 0, Zenv=50% (solid line) differ by less 
than 4%. Therefore, for such levels of Zonv, the models are 
less sensitive to the treatment of the Z-element in the enve- 
lope. As mentioned previously, the sequence based on the 
ANEOS EOS (short-dash line) yields a significantly faster 
cooling sequence, with a ~ 6% — 10% smaller radius at 1 
Gyr, compared with all other sequences. This large uncer- 
tainty on the EOS will be the major culprit for preventing 
accurate determinations of the exoplanet internal composi- 
tion from their observed radius. 

We have also explored the effects of the heavy element 
distribution for smaller total enrichments, Z < 15%. In this 
case, we find a < 2% global effect on the radius at a given 
age, depending whether the heavy elements are all located 
in the core or are distributed throughout the entire planet, 
using either an effective He abundance or the AVL with 
the SESAME EOS for the thermodynamics of the heavy 
component. Models based on the ANEOS EOS (with heavy 
elements distributed over the entire planet) still predict the 
smallest radii at a given time, with a maximum ^ 10% 
effect compared with the other sequences. 

Figure [5] portrays the internal T{P) and p{P) profiles 
for the various models of our 20 M® planet at an age of 4 
Gyr. 



The first conclusions to be drawn from these tests for a 
Neptune-mass planet are the following: 

(i) for a core mass less than 50% of the planet's mass, a 
variation of the core composition from pure water to pure 
iron, the maximum expected difference in mean densities, 
yields a difference on the radius of less than ~ 10% after 1 
Gyr. Yet, such a difference is accessible to the observational 
determination of some transit planet radii, if (when) other 
sources of uncertainties, in particular on the heavy material 
EOS, were (will be) under control. 

(ii) for a metal- fraction in the envelope Zenv < 20%, 
the EOS of the Z-material can be approximated by using a 
corresponding Ycquiv effective helium fraction in the SCVH 
EOS. This approximation become more dubious, and even 
wrong, above this limit. 

(iii) if Z <^15%, either globally or in the envelope, the 
different treatments of heavy elements yield a relative vari- 
ation of the radius of ^ 2%, except when using the ANEOS 
EOS, which yields a ~ 10% difference. Within this limit 
for Z, and given all the other uncertainties in planet mod- 
elling, the impact of heavy elements on the evolution of 
Neptune-mass planets can be mimicked reasonably well by 
considering that the heavy elements are all located in the 
core. 

(iv) For larger heavy element enrichments (Z>15%), 
the distribution (everything in the core versus uniform dis- 
tribution) and the treatment of the heavy element thermo- 
dynamic properties ( Yoquiv or AVL) can significantly affect 
the cooling and thus radius determination for a given age 
(more than 10% difference after 1 Gyr). These uncertain- 
ties, unfortunately, hamper an accurate determination of 
the detailed composition of the heavy element material in 
the planet's interior. 

4.2. Jovian-mass planets 

In this section, we extend the analysis done in the previ- 
ous section to a template 1 Mj planet. We first analyse 
the effect of core composition for a core mass of 159 Mg, 
corresponding to a total heavy element enrichment Z=50% 
(see Fig. |B^). We find slightly larger effects than for the 
Neptune-mass case, with a 7% (15%) difference in radius 
between the pure water and the pure rock (iron) core cases. 
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Fig. 3. Effect of the composition and internal distribution 
of heavy elements on the radius evolution for a planet of 
mass Mp = 20 with a total heavy element mass fraction 
Z = {Mz/Mp)=50%. (a): Models with Mcoro = lOMe 
of water (solid line), rock (dashed-line) and iron (dotted 
line), and Zcnv=0. (b): Models with no core and heavy el- 
ements distributed over the entire planet. Dash-dotted line: 
mixture of H/He + Z described by the SCVH EOS with 
Y equiv— 0-275 -|- 0.5 — 0.775. Long-dashed line: mixture of 
H/He + water (SESAME EOS) using the additive volume 
law. Short-dashed line: mixture of H/He + water (ANEOS 
EOS) using the additive volume law. For comparison, the 
model with Mcore = 10 M^ of water and Zcnv=0 is shown 
by the solid line, (c): More realistic models with a 8 M^ 
core of water (ANEOS EOS) and a Ze„v = 16.6 % heavy ele- 
ment enrichment in the envelope. Dash-dotted line: mixture 
of H/He + Zenv described by the SCVH EOS with Yoquiv= 
0.275 + 0.166 = 0.441. Long-dashed line: mixture of H/Hc 
-I- water (SESAME EOS) using the additive volume law. 
Short-dashed line: mixture of H/He -I- water (ANEOS EOS) 
using the additive volume law. Also displayed is the model 
with Mcoro = 10 M® of water and Zonv=0 (solid line). 

respectively. Whether such massive cores can indeed form 
for Jovian-mass planets will be considered in §7.1. As for 
the Neptune case, using either the ANEOS or the SESAME 
EOS for the same material in that case yields similar cool- 
ing sequences, with < 3% differences on R for water and 
< 2% for rock. For these planets, the iron core does not un- 
dergo a phase transition and iron always remains liquid, ac- 
cording to the ANEOS EOS. Indeed, central temperatures 
and densities are significantly larger than for Neptune-mass 
planets, with Tc ~ 10^ K and pc ^ 157 g.cm"-^ at 1 Gyr for 
a 1 Mj planet with an iron core (see Table 2). 




0.2 0.4 0.6 0.8 1 

m/Mp 



Fig. 4. Effect of the EOS on the inner profile as a func- 
tion of fractional mass for a planet of 20 M® with a 
mass fraction Z=50% of heavy elements distributed over 
the entire planet. In all panels, the curves correspond to: 
a mixture of H/He + Z described by: the SCVH EOS 
with Yequiv=0.275 + 0.5=0.775 (solid line); a mixture of 
H/He + water (SESAME EOS) using the additive vol- 
ume law (dash-dotted line); a mixture of H/He -I- water 
(ANEOS EOS) using the additive volume law (dashed line). 
Upper panel: Inner profile of the local gravothermal en- 
ergy, —TdS/dt (in erg/g/s), at the beginning of the evo- 
lution. Middle panel: Pressure profile (in dyne/cm^) at 1 
Gyr. Lower panel: Specific entropy profile (in erg/g/K) at 
1 Gyr. 



We explore the effects of the heavy element distribu- 
tion with Z = 50% (Fig. O)) and Z = 20% (Fig. W)- The 
latter case is comparable to the expected level of enrich- 
ment in Jupiter and Saturn (Saumon & Guillot 2004). In 
the Z — 50% case, both the various thermodynamic treat- 
ments and localizations of the heavy elements yield signif- 
icantly different cooling behaviors (see Fig. Wp)- The se- 
quence based on the Yoquiv approximation, notably, differs 
from the other ones, as expected from our previous study 
for Neptune-mass planets. At 1 Gyr, the radii obtained with 
models based on the AVL with ANEOS and SESAME, re- 
spectively, differ by ~ 6%. The different distributions of 
heavy elements within the planet have an even larger im- 
pact, with up to 12% difference in R between the case with 
Mcore=159 M0, Zonv=0 and the case with Z distributed 
throughout the entire planet with ANEOS. Interestingly 
enough, for the present Jovian conditions, models based 
on the ANEOS or SESAME EOS, when metals are mixed 
throughout the entire planet, yield smaller differences than 
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Table 2. Radius and central thermodynamic properties of a planet of mass Mp = lAfj at 1 Gyr for two heavy element 
mass fractions Z = {Mz /Mp)=50% and Z = 20%. The labels are the same as in Table [TJ As mentioned at the end of 
§4.2, for Z=50%, the solution with Mcore=0 and Z = Zenv is similar to a solution with Mcore ~ 10 and the rest of 
the heavy material distributed in the envelope. 
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Fig. 5. Internal temperature (K) and density (g cm^'^) pro- 
files for a 20 planet, irradiated by a Sun at 0.045 AU, 
at an age of 4 Gyr. Solid line: Mcore = 0.5 Mp, Zonv=0; 
long-dash line: M^ore = 0, Z = 0.5, AVL SESAME EOS; 
short-dash line: M^ore ^ 0, Z ^ 0.5, AVL ANEOS EOS. 



for the Neptune-mass planet case, after about 1 Gyr, i.e. 
at high pressure and moderately high temperature. Indeed, 
the inner pressure and temperature conditions are very dif- 
ferent between Neptune-mass and Jupiter-mass planets (see 
Figures [5] and El) . For the former ones, P ranges between 



10"'^° and 10^^ dyne.cm"^ over more than 99% of the total 
mass (see Fig.S]) while T ranges from ^ 5000 to lO"' K. For 
the latter ones, the typical domains are P = 10^^ — lO^** 
dyne.cm-2 and T = 10"* - 10^ K. As mentioned in §3.2, dif- 
ferences in the P and T dependence of the entropy between 
the different EOS are more pronounced under the interior 
conditions of Neptune-mass planets than for Jupitcr-mass 
planets (see Fig. 

If, instead of the extreme Mcore = or Zcnv = cases, 
we take a more realistic model with a 10 core, com- 
parable to what is expected in Jupiter or Saturn, and we 
distribute the rest of heavy elements homogeneously in the 
H/He rich envelope, we find essentially the same evolution 
as when the heavy elements are distributed throughout the 
whole planet, with no core. This suggests that for mas- 
sive, metal-rich planets, the evolution is better described 
by models which assume that all metals are distributed 
over the entire planet, since this yields results similar to 
the ones obtained with a more realistic distribution, than 
by models which assume that all heavy elements are in the 
core, with a metal-free, Z — envelope. 

For a more moderate heavy element enrichment, Z = 
20% (see Fig. 15};) , the treatment of heavy elements in the 
entire planet, based on an Ycquiv or on the AVL with 
ANEOS or SESAME, is found to be less consequential, for 
the present Jupiter-mass planet case. The different meth- 
ods to describe the EOS yield less than 2% differences on 
i? at a given age. The effect of the distribution of heavy 
elements (all in the core versus all distributed over the en- 
tire planet) is slightly more consequential, with up to 4% 
difference on R. 

Figure [7] portrays the internal T{P) and p{P) profiles 
for the various models of our 1 Mj planet at an age of 4 
Gyr. 

The conclusions derived from these tests for a Jupiter- 
mass planet can be summarized as follows: 

(i) for a core mass less than 50% of the planet's mass, a 
variation of the core composition between pure water and 
pure rock (iron) yields an effect on the radius of < 7% 
(15%). 
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Fig. 6. Effect of tlie composition and the distribution of 
heavy elements on the radius evolution of a planet of 1 A/j 
(318 Mq). (a): Total heavy element enrichment Z — 50%. 
Models with Mcore = 159 of water (solid hne), rock 
(dashed-hne) and iron (dotted hne), and Zenv=0. (b): 
Total heavy element enrichment Z — 50%. The solid line 
corresponds to a model with Mcore = 159 Mq of water 
and Zcnv=0. The other curves correspond to models with 
no core and heavy elements distributed over the entire 
planet: mixture of H/He + Z described by the SCVH EOS 
with Yequiv=0.275 + 0.5=0.775 (dash-dotted line); mix- 
ture of H/He water (SESAME EOS) using the additive 
volume law (long-dashed line); mixture of H/He -I- water 
(ANEOS EOS) using the additive volume law (short-dashed 
line), (c): Total heavy element enrichment Z — 20%. Solid 
line: model with Mcore — 63.6 M® of water and Zenv=0; 
dash-dotted-line: no core and mixture of H/He -I- Z de- 
scribed by the SCVH EOS with Ycquiv=0.275 0.2 = 
0.475; long-dashed line: no core and mixture of H/He -|- 
water (SESAME EOS) with the AVL and Zo„v=0.20; short- 
dashed line: no core and mixture of H/He + water (ANEOS 
EOS) using the AVL and Zc„v=0.20. 



(ii) for a global metal-enrichment Z <20%, the treat- 
ment and distribution of heavy elements affect only mod- 
estly the radius predictions (less than 4% effect on R at a 
given age). 

(iii) For significant heavy element enrichment (Z ^ 
50%), the distribution of heavy elements (everything in the 
core versus uniform distribution) has a non negligible effect 
on the radius predictions, up to 12% difference at a given 
age. For such high metal-enrichments, distributing all the 
heavy material throughout the entire planet, with no core, 
is found to yield similar results as a more realistic model 



Fig. 7. Internal temperature and density profiles for a 1 



A/j planet at an age of 4 Gyr. Solid line: Mc. 
long-dash line: M, 
short-dash line 



0.5 M, 



core - 0, Z = 0.5, AVL SESAME EOS; 
Mcore = 0, Z = 0.5, AVL ANEOS EOS. 



with about a 10 M^ core and the rest of the material in the 
envelope, while a model with a metal-free envelope and all 
heavy elements in the core differs substantially from these 
sequences. 

4.3. Energy release from the core 

As mentioned in §2.2, the thermodynamics of the heavy 
material component in planet cores is often described either 
by a zero-temperature EOS or under the assumption of a 
uniform, low temperature (Fortney et al. 2007; Seager et 
al. 2007), or a uniform density (Bodenheimer et al. 2003). 
These models ignore the gravothermal energy contribution 
of the core, (— T^)core to the cooling of the planet. As 
an other extreme. Burrows et al. (2007) arbitrarily assume 
that the specific entropy of the core material is the same 
as that of the H/He envelope, therefore overestimating the 
heat release of the dense core. 

According to the ANEOS or SESAME EOS, the typical 
heat capacity of water in the cores of the presently studied 
Neptune-mass and Jovian-mass planets is Cv ~ 3 x 10^- 
5x10^ erg g~^ K~^, for typical temperatures T ^ 5000- 
5x10^ K and pressures P ~ 10^^ — 10^^ dyne cm~^ (see 
Figs. [5] and [7|). This is about 1/3 the specific heat of 
H/He in the envelope. This corresponds to variations of 
the specific heat from the high-temperature nearly ideal gas 
limit, ^ 6TZ/A (where A is the species mean atomic weight 
and TZ is the perfect gas constant), to regimes which in- 
clude potential-energy contributions associated with trans- 
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lational and librational modes (3TZ/A maximum for three 
translational and three hbrational modes per molecule). If 
the species enters the solid phase, the specific heat decreases 
rapidly (Debye regime) and eventually vanishes. The spe- 
cific heat of rock or iron is even smaller (larger atomic mass) 
and the contribution of cores made of such materials to the 
cooling history of planets is negligible. For planets in the 
mass range 20 - 1 Mj, as explored in the previous sec- 
tion, with core masses less than 50% of the planet's mass, 
and with cores made up of water, the release of gravitational 
energy (P^) and of internal energy (^) never exceeds 40 
% of the total released energy, T{^). For Neptune-mass 
planets, the dominant contribution from the core is due to 
its contraction, i.e the work, during the first 0.5 - 1 

Gyr of evolution, and by the release of its internal energy at 
older ages. For jovian-mass planets, the contraction of the 
core provides the dominant contribution to its gravother- 
mal energy during the entire planet evolution. 

The assumption of zero-T or uniform low temperature 
affects the structure of the cores, because of the temper- 
ature dependence of the P{p) relations predicted by cur- 
rent EOS (see §3.2 and Fig. [T|). But, more importantly, 
it affects the evolution of the planet, because it implies a 
negligible contribution of the core to the total gravother- 
mal energy of the planet. Indeed, the release of thermal 
energy from the core, (^)corc, is forced to be zero, and 
the small compressibility of water at low T drastically un- 
derestimates the true contraction of the core during the 
planet evolution, and thus the release of gravitational en- 
ergy. We have analysed this assumption by imposing a con- 
stant and uniform temperature of 300 K in the water core 
of the planet models analysed in the previous section, us- 
ing ANEOS and SESAME. For the jovian-mass planet (1 
Mj) with Mcore=159 M^, the effect on the radius is neg- 
ligible (less than 2%) after 1 Gyr. On the Neptune-mass 
planet (20 M^), the effect is larger, the largest effect being 
found with the SESAME EOS for the model with a 10 M® 
core, that yields a 6% effect on R after 1 Gyr. The largest 
temperature variations of the EOS are indeed found at low 
pressure (P < 1 Mbar), as illustrated in Fig. [TJ affecting 
more importantly light planets. 

The heat transport in the core is due to convection or 
conduction (by electrons or phonons), depending on the 
core material and the age of the planet. Conduction may 
dominate after a few Gyr of evolution, as the core becomes 
cool and dense enough for the thermal conductivity due to 
degenerate electrons to become large enough, Kc oc p^/'^, 
and the conductive flux to dominate the convective one, 
J'cond = —Kc^T > J-conv It starts to dominate earlier 
for rock and iron cores, compared to water cases. These 
results, however, are hampered by the uncertainties in the 
conductive opacities calculated with the ANEOS EOS or 
with Potekhin (1999) for such materials. We have checked 
the effect of such an uncertainty on the cooling history in 
the case of the 20 M^ planet with Mcorc=10 of water 
with the SESAME EOS, as this sequence provides a case 
with the largest energetic contribution from the core, for en- 
richments Z < 50 %. Over the entire evolution, heat is pre- 
dicted to be transported by convection within the core and 
the temperature gradient is given by the adiabatic gradient. 
We have arbitrarily decreased the conductive opacities so 
that conduction now dominates over convection. The core 
thus becomes isothermal. The effect remains small on the 



radius evolution (maximum 3% on i? at a given age com- 
pared to the convective case). Therefore, the uncertainty 
in the heat transport efficiency of the core has a smaller 
impact than neglecting the temperature dependence of the 
core material and its global energetic contribution. 

To conclude this section, the present calculations show 
that neglecting the thermal and gravitational contributions 
from the core, by assuming zero-T or low uniform T, in cur- 
rent planet modelling, leads to a maximum effect of ~ 6% 
on the radius, after 1 Gyr of evolution, for cores less massive 
than 50% of the planet's mass. Whether such a variation 
can be considered negligible or not depends on the accuracy 
of the data the models need to be confronted to in order 
to infer the planet's internal composition, once the EOS of 
heavy elements will be determine with enough accuracy. As 
will be shown in the next section, the conclusion is dras- 
tically different for larger metal contents, for Neptune-size 
planets. 

5. Extreme metal-enrichment: the two specific 
cases of HD149026b and GJ436b 

In this section, we analyse extreme cases of heavy element 
enrichments (Z > 50%), focusing on two recently discov- 
ered transiting planets, the Saturn-mass planet IID149026b 
and the Neptune-mass planet GJ436b. Our goal is to ex- 
amine the properties of these two transit planets, in the 
light of the analysis conducted in the previous section con- 
cerning the main assumptions and uncertainties of current 
planetary models. 

5.1. HD149026b 

The discovery of the Saturn-mass planet HD149026b (Sato 
et al. 2005) revealed an unexpectedly dense planet, with 
a mass Mp = 0.36 ± 0.04 Mj and a radius Rp = 0.725 ± 
0.05 i.e. a mean density p = 1.17 ± 0.35 g.cm~^. For 
comparison, Saturn has a mass of 0.3 Mj but a radius of 
0.8 i.e. p — 0.66 g.crn"^. The transit planet is orbiting 
a G-type star at an orbital distance a = 0.042 AU, about 
230 times closer to its Sun than Saturn. The age of the 
system, ~ (2 ± 0.8) Gyr, is also younger than our Solar 
System. It is important to understand the nature and the 
origin of this puzzling planet, and in particular to know 
whether current planet formation scenarios, in particular 
the core accretion model, can explain it. This requires a 
knowledge of its structure and composition, which can only 
be inferred from theoretical models. Several authors have 
tried to infer the inner structure of this planet under the 
usual assumptions described in §2. The models assume that 
almost all heavy elements are in the core, and the H/Hc 
envelope is either free of metals (Burrows et al. 2007), or is 
moderately enriched, with Z = Zq (Ikoma et al. 2006) or 
Z = 0.045 (Fortney et al. 2007). Ikoma et al. (2006) have 
also investigated a case with Zonv = 0.37. Fortney et al. 
(2007) use a zero-temperature EOS for the core and thus 
ignore its heat content contribution. They find negligible 
effect when using a nonzero temperature EOS for the ice 
mixture given by Hubbard & Marley (1989). This latter 
EOS provides pressure-density relations appropriate for the 
description of warm adiabatic mixtures (T ^ 10"*) K, but 
do not explicitly account for the temperature dependence. 
Ikoma et al. (2006) also use this EOS and assume a uniform 
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Fig. 8. Evolution of a planet with the characteristics of 
HD149206b (mp=116 ~ 0.36Mj) with different heavy 
element distributions. Solid line: model with Mcorc — 
80 Me of water (ANEOS EOS) and Zc„v=0. The other 
curves correspond to models with heavy elements dis- 
tributed over the entire planet. Dash-dot line: mixture of 
H/He + Z described by the SCVH EOS with an effective 
helium abundance, Ycquiv=0.275-f 0.69=0.965, which cor- 
responds to Z = 0.69, i.e 80 of heavy elements. Long- 
dash line: mixture of H/He -I- water (SESAME EOS) using 
the AVL with Z=69%. Short-dash hue: mixture of H/He + 
water (ANEOS EOS) using the AVL with Z=69%. Dotted 
hne: Mixture of H/He + water (ANEOS EOS) using the 
AVL with Z=51%, i.e 60 M^ of heavy elements. 

temperature in the core and a specific heat Cv ~ 10^ erg 
g~^ to account for the core heat release. According to 
all these models, the total mass of heavy elements in the 
planet hes in the range ~ 40-90 M®, i.e Z ~ 35% - 80%. 

We will test the impact of these main assumptions done 
in current structure models of HD149026b, namely the dis- 
tribution of all heavy elements in a core and the use of 
zero-T EOS. Since the present study focusses on the sen- 
sitivity of the structure and evolution for a given setup of 
atmosphere models, we use, as outer boundary condition to 
the inner structure, atmosphere models with solar compo- 
sition, Zq, irradiated by a G-star at 0.04 AU, even though 
the results are expected to change to some extent with the 
composition of the atmosphere (see e.g Burrows et al. 2007). 
This issue will be explored in a forthcoming paper. We ob- 
tain a good match of the planet's radius with a model with 
a water core of 80 Mq and a metal-free envelope, Zcnv — 0. 
Figure [8] displays this model (solid line), calculated with the 
ANEOS EOS in the core. The use of the SESAME EOS for 
water in the core yields a 3% smaller radius, still providing 
a good match to the observed value. Figure [8] also shows 
the impact of the heavy element distribution, with models 
where the 80 Mq of heavy elements are distributed over the 
entire planet. As anticipated from the studies conducted in 
§3.2, the use of an Yoquw in the SCVH EOS to handle the 
heavy element contribution yields, for such a high metal 
fraction, an evolution that differs drastically from the other 



ones, even though, coincidentally, it gives a good solution at 
the age of the system. As expected from the tests performed 
in §4, models with 80 Mq of heavy elements mixed with 
H/He throughout the entire planet, using the AVL with 
the SESAME (long-dash hne) or the ANEOS (short-dash 
line) EOS for water, yield denser structures. The latter one 
yields a radius ^ 25% smaller than the observed value, an 
effect larger than changing the core composition from water 
to rock (~10% effect on R), for the same core mass of 80 
Me, and similar to the one obtained when changing from a 
pure water to a pure iron core, a rather unlikely solution. 

Instead of the two aforedescribed extreme heavy ele- 
ment distributions, we have also examined a more realistic 
model with a 20 M® core, similar to what is predicted for 
Saturn, and the remaining 60 M® mixed within the H/He 
envelope, using the AVL both with SESAME and with 
ANEOS. As anticipated from the tests performed in §4.2, 
this yields similar results than assuming that all heavy ele- 
ments are distributed throughout the entire planet. Models 
with such metal enrichment and distributions thus seem to 
be excluded by the observations. An alternative solution for 
HD149026b, however, is obtained with a 60 Mq of heavy 
elements, instead of 80, mixed with H/He throughout the 
entire planet, described with the ANEOS EOS for water 
and the AVL formalism (dotted line in Fig. [8]), whereas 
when using the SESAME EOS for water, 70 M© are re- 
quired. As mentioned above and as shown in our previ- 
ous studies, such models are equivalent to models with a 
small core mass and a significantly enriched envelope. This 
shows that for a given heavy element material, water in the 
present case, the effect of modifying its internal distribu- 
tion (everything in the core or a fraction redistributed in 
the H/He envelope) has by itself a large impact, yielding 
an uncertainty on the amount of heavy material required to 
reproduce the observed radius of 80-60=20 M©. As shown 
in the previous sections, assuming all heavy elements to be 
in the core yields structures less dense than when these el- 
ements are mixed throughout the envelope, so that these 
models require a larger amount of heavy material to match 
a given radius. 

Finally, as done in §4.3, we have tested the effect of 
the temperature dependence of the EOS in the core by as- 
suming a uniform temperature of 300 K in a model with a 
80 Me water core, using ANEOS and SESAME EOS. The 
largest effect is found with this latter EOS, the model with 
a uniform low temperature yielding a radius 3 % smaller 
than the "hot" core case at 2 Gyr. 

5.2. GJ436b 

The first Neptune-mass transiting planet has been discov- 
ered recently (Gillon et al. 2007a), with a mass Mp = 
22.6 ±1.9 Me and a radius Rp = 25200 ±2200 km = 0.352 
±0.03 The planet is orbiting an M-type star of ~ 0.44 
Mq at an orbital distance a ^ 0.028 AU. According to the 
models of Fortney et al. (2007), Gillon et al. (2007a) sug- 
gest that the planet is composed predominantly of ice with 
a thin H/He envelope of less than 10% in mass. Based on 
Spitzer observations, Gillon et al. (2007b) and Deming et 
al. (2007) determine a shghtly larger radius, with i?p=0.386 
±0.016i?j. We adopt this value in the following and we as- 
sume an age for the system of 1-5 Gyr, as it is essentially 
unconstrained by the observations. Models available in the 
literature to determine the inner structure of this planet 
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use temperature-independent EOS to describe the core and 
thus ignore its thermal contribution (Fortney et al. 2007; 
Seager et al. 2007; Adams et al. 2007). We will test this 
assumption. 

We have calculated models characteristic of GJ436b. 
We use solar metallicity atmosphere models and, given the 
low luminosity of the parent star, we neglect presently the 
irradiation effects. Recent determinations of GJ436b irra- 
diation induced temperature (i.e. brightness temperature) 
indeed suggest that the evolution is not likely to be signifi- 
cantly altered by irradiation (Demory et al. 2007). A model 
with a core of 21 made of water and a metal- free en- 
velope, Zonv=0, provides a good match to the observed 
radius, as illustrated in Fig.[9](sohd and dashed lines). The 
ANEOS EOS for water (solid line) yields a slightly larger 
(~ 4%) radius than the SESAME EOS (dashed hne). A 
good match is also obtained with a rocky (dunite) core of 
mass 19.5 M^. Given the high mass fraction of heavy ma- 
terial of this planet, more than 80%, the freedom to vary 
its distribution is limited. Assuming a core of 21 M0 with 
Zenv=0 or a slightly smaller core with some heavy element 
enrichment in the envelope, a more realistic solution, e.g. 
a model with Mcore = 20 M0 and Zonv=0.38, yield less 
than 8% variations on the radius, still within the observed 
error bars. According to the tests performed in §4.1 on a 
20 M0 planet, larger values of Zenv ( Zenv ^ 50%), and 
thus smaller cores, yield larger effects on the radius. Note 
that, according to present planet formation models (Alibert 
et al. 2005a; Mordasini et al. 2008), such high fractions of 
Z-material enrichment in envelopes of Neptune-mass plan- 
ets are not excluded (Baraffe et al. 2006, figure 1) since, as 
mentioned in §2.1, the fate of accreted planetesimals during 
the gas accretion phase depends on the envelope mass. 

Figure [9] also displays the evolution of the aforemen- 
tioned model with the assumption of a uniform cold tem- 
perature (300K) in the core. This yields a moderate effect 
on the evolution, about 4% on the radius at 1-5 Gyr with 
the ANEOS EOS and - 6% effect with SESAME. In the 
latter case, however, the model lies outside the observa- 
tional error bars (see dotted line in Fig. [9]). As seen in Fig. 
[TUl for planets with such a large fraction of heavy mate- 
rial, and conversely with such a modest gaseous H/He frac- 
tion (^ 10%), the evolution of the planet is largely domi- 
nated by the core contribution. Neglecting the temperature- 
dependence of the core, and thus its contribution to the 
gravothermal energy of the planet, yields an incorrect cool- 
ing sequence. Even though, for such low-entropy material, 
the global gravothermal energy, TdS/dt, remains small, a 
proper calculation should take into consideration the core 
contribution. Note also that even a modest H/He fraction 
affects the radius determination, as illustrated in Table [3] 
by comparing the present results with the value of the ra- 
dius corresponding to a pure 22.6 M^ icy planet with no 
gas envelope, as derived from the fitting formulae of Seager 
et al. (2007). 

Importantly enough, taking into account the thermal 
and gravitational energetic contributions from the core be- 
comes even more crucial if the irradiation effects from 
the parent star are important. Fig. [5] shows the evolu- 
tion of the same planet model, with Mcorc = 21 of 
water (SESAME EOS), but with an incident stellar flux 
^inc = sl"^)^-^*' where is the flux from the parent 
star, six times larger than for GJ436a. This corresponds to 
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Fig. 9. Evolution of a planet characteristic of GJ436b 
(A/p=22.6 M^,) with different heavy element distributions. 
Models with no irradiation: Solid line: model with Mcorc — 
21 Me of water (ANEOS EOS) and Zo„v=0; short-dash 
line: model with Mcorc = 21 M© of water (SESAME EOS) 



and Z„ 



=0; dotted line: test model with Mr. 



21Mfl 



of water (SESAME EOS) and Zenv=0, assuming an uni- 
form temperature of 300 K in the core. Effect of irradia- 
tion, with Fine = 6 X Fine for GJ436b (see text): long-dash 
line: same as short-dash line with irradiated atmosphere 
models; dash-dot line: same as dotted line with irradiated 
atmosphere models 
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Fig. 10. Various contributions to the gravothermal energy. 
Eg = —TdS/dt, normalized to the total value, for a planet 
characteristic of GJ436b (Mp=22.6 M©), with a 21 M© 
core of water with the SESAME EOS. Sofid line: global 
contribution from the core; dash-dot line: contraction work 
contribution from the core; long-dash line: thermal contri- 
bution from the core; short-dash line: global contribution 
from the H/He 1.6 M© envelope. 
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Table 3. Radius of a 22.6 planet at 2 Gyr with dif- 
ferent water core masses, EOS and levels of irradiation, 
characterised by a stellar incident flux, Fi„c- Two cases are 
considered: no irradiation (Finc=0) and Finc=6 x Fine of 
GJ436b. Comparison is done with the pure water (ice) case 
of Seager et al. (2007). Results based on the assumption of 
a uniform core temperature of T=300K are also given. 
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a parent star about 50% hotter. In that case, neglecting 
the temperature dependence of the EOS in the core and its 
contribution to the planet's cooling yields a ~ 11% (7%) 
smaller radius with SESAME (ANEOS) EOS (dash-dot vs 
long-dash curves). This stems from the larger planet in- 
terior temperature and entropy in the irradiated sequence 
compared with the non-irradiated one. In the irradiated 
case, the core temperature ranges from 5000 to ~ 2 x 10^ K 
(from the bottom of the H/He envelope to the center), to be 
compared with 3000 to ^ 10"* K in the non-irradiated case, 
characteristic of the temperatures expected in the rocky/icy 
part of Neptune and Uranus (Guillot 2005) . 

6. Evolution of super Jupiter planets: Hat-P-2b 
and deuterium burning planets 

The final part of our study is devoted to "super- Jupiter" 
extra-solar planets, with masses Alp ^ 1 Mj. It is mo- 
tivated by the growing number of discoveries of massive 
extra-solar planets, in the mass regime overlapping the one 
of low-mass brown dwarfs, issued from a different forma- 
tion mechanism. These discoveries feed the heated debate 
concerning the definition of a planet and the possibility 
to distinguish planets from brown dwarfs of similar mass. 
In this context, one of the most remarkable discoveries is 
the transiting super-Jupiter object HAT-P-2b (also named 
HD147506b), with a mass Mp = 9.04 Mj and a radius Rp 
= 0.982 Rj (Bakos et al. 2007). Loeillet et al. (2007) reanal- 
ysed the orbital parameters of the system and find similar 
values, Mp = 8.64 Mj, Rp = 0.952 Rj. In that case, the 
mass-radius relationship offers a unique information to infer 
the gross composition of the object and to determine its real 
nature, low-mass gaseous brown dwarf or very metal-rich 
massive planet. Bakos et al. (2007) suggest that the mean 
density of this planet is only marginally consistent with 
model predictions for an object composed predominantly 
of H and He, and requires the presence of a large core, with 
Mcore > 100 Mq. Here, we calculate more thoroughly the in- 
ternal structures consistent with the radius determination, 
along the lines described in the previous sections. We find 
that a total amount of - 300-600 (200-500) of a water 
(rock) component is required to explain the radius at the 
present age, as illustrated in Fig. [TlJ The rather large pre- 
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Fig. 11. Evolution of a super- Jupiter planet with the char- 
acteristics of HAT-P-2b (Mp= 9 Mj) for different distribu- 
tions of heavy elements. Irradiation effects from the par- 
ent star are taken into account, adopting a time- average 
orbital distance a — 0.077 AU, for the proper eccentricity 
e=0.52 and semimajor axis of the relative orbit arei=0.0677 
AU. SoHd line: Mcoro = 600 M® of water (ANEOS EOS) 
and Zcnv=0. Long-dashed line: Mcore — 350 of wa- 
ter (ANEOS EOS) and Ze„v=0. Dash-dotted hne: 600 
of heavy elements with 100 M® in the core (water, with 
ANEOS EOS) and 500 M® in the envelope, i.e. Z^nv = 
0.18, mimicked by an Yoquiv=0.275 -I- 0.18. Short-dashed 
line: H/He brown dwarf with Z = Zq. 

dieted range of heavy material enrichment stems from the 
large observational error bars. As seen in the figure, a brown 
dwarf (Baraffe et al. 2003), i.e. a gaseous H/He object with 
solar metallicity, is predicted to have at this age a radius 
marginally consistent at the 2a limit with the observations 
(short-dashed line) and thus can be excluded at the ^ 95% 
confidence level, if the present radius observational determi- 
nation is confirmed. This shows that planets can exist with 
masses up to ~ 10 A/j, well above the opacity-limit for frag- 
mentation, m ~ 3-5 Mj (jWhitworth fc Stamatellos 2006p . 
the expected minimum mass for brown dwarf formation. 
These two distinct astrophysical populations should then 
overlap over a substantial mass domain. 

In this context, it is interesting to explore the fate of 
planets massive enough to ignite deuterium-fusion in their 
central parts, i.e. with Mp>12Mj (Saumon et al. 1996, 
Chabrier et al. 2000). Indeed, recent calculations of plane- 
tary population synthesis based on improved core-accretion 
models of planet formation (Alibert et al. 2005a) and a wide 
variety of initial conditions, predict the formation of super- 
massive planets, up to ^ 20-25 Mj, with rocky/icy core 
masses up to several 100 M® (Ida & Lin 2004; Mordasini 
et al. 2007, 2008). If planets with a massive core can 
form above the aforementioned deuterium burning mini- 
mum mass, a key question is to determine whether or not 
the presence of the core can prevent deuterium burning to 
occur in the deepest layers of the H/He envelope. Guided 
by the results of Mordasini et al. (2008), we have consid- 
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ered a 25 Mj planet with a 100 core. Independently 
of the composition of the core material (water or rock), 
deuterium-fusion ignition does occur in the layers above the 
core and deuterium is completely depleted in the convec- 
tive H/He envelope after ~ 10 Myr. The same conclusion 
holds for a core mass of several 100 M®. These results high- 
light the utter confusion provided by a definition of a planet 
based on the deuterium-burning limit. 

7. Planet evolutionary models and mass-radius 
relationships 

7.1. The reality test 

Some of the internal structures determined in the previous 
sections in order to match the observed radii and inferred 
mean densities of transit planets are rather unusual, and 
the possibility to form such compositions must be exam- 
ined in the context of our current understanding of planet 
formation. According to current models of planet forma- 
tion, which include migration (Alibert 2005a, 2006;), up to 
~ 30% of the heavy material contained in the protoplane- 
tary disk can be incorporated into forming giant planets 
(Mordasini et al. 2007). The maximum mass for a sta- 
ble protoplanetary disk is Md <0.1 Mi,, so that, for a 1 
Mq parent star of solar composition, Z ~ 2%, as much as 
Mz « 0.3 X 0.1 X 0.02 X (3.3 x 10^) w 200 of heavy 
material can be accreted onto the planet. Present planet 
formation models (Mordasini et al. 2007) reach about this 
limit for very massive planets, for a 1 Mq parent star. 
Therefore, in principle, according to these calculations, a 
heavy element mass fraction > 50% can not be excluded, 
even for Jovian-type planets. This requires, however, accre- 
tion rates from the planet's feeding zone significantly larger 
than the values typical of the early "runaway growth" ac- 
cretion phase, ~ 10~^ yr~^ (Ward 1996). It also im- 
plies migration of the planet's embryo, or some substantial 
orbit eccentricity. Indeed, in the absence of migration or 
eccentricity, tidal interactions between the planet and the 
disk are supposed to lead to the opening of a gap once 
the planet has reached about a Saturn total (gas-|-solid) 
mass, ~ 100 M^, for the minimum mass solar nebulae con- 
ditions, after which planetesimal accretion decreases dra- 
matically (Lin & Papaloizou 1986). Note that, if the plan- 
ets formed originally at large orbital distances and migrated 
inwards, they are expected to have a significant content of 
heavy material, given the larger available mass reservoir. 
Furthermore, short-period planets are expected to have a 
larger fraction of heavy material than planets located fur- 
ther away for two reasons. First of all, more impacts from 
low eccentricity orbit planetesimals are expected. Second 
of all, the closer the orbit of the planet, the larger its or- 
bital speed, (GM*/a)^/'^, compared to its escape velocity, 
(2GMj,/i?^)^/^, making planetesimal ejection less efficient 
(Guillot 2005). 

A second concern is the possibility to have the type of 
heavy element distributions examined in the present paper. 
As mentioned earlier, the fact that gas accretion is triggered 
once the core mass has reached about 6 to 10 Mq seems 
to be a rather robust result (Mizuno 1980, Stevenson 1982, 
Pollack et al. 1996, Alibert et al. 2005a, Rafikov 2006). For 
larger enrichments, the rest of heavy elements should thus 
be mixed with the H/He envelope. As shown by Stevenson 
(1982), the maximum amount of heavy material (compared 



with the H/He medium) which can be redistributed uni- 
formly by convection throughout the planet from an ini- 
tially stably stratified configuration is of the order of the 
planet's mass. So in principle, a planet with no or small core 
but all or most of the heavy material being redistributed 
throughout the gaseous envelope is possible. As shown in 
our study, for Jovian type planets (see §4.2), for a global 
metal fraction Z = 50%, if we assume a core of 10 Mq, 
comparable to what is expected in Jupiter or Saturn, and 
we distribute the rest in the envelope, we find essentially 
the same evolution as when the heavy elements are dis- 
tributed throughout the whole planet, with no core. This 
supposes, of course, that large-scale convection remains efh- 
cient enough to redistribute homogeneously the heavy ma- 
terial in the H/He envelope (Chabrier & Baraffe 2007). Note 
that this result also holds for lower mass planets with simi- 
lar enrichment, Z = 50%, if the core mass remains < 10% 
Mp. 

Therefore, the enrichment in heavy material and the 
internal compositions explored in the present calculations 
have at least some reasonable theoretical foundation and 
can not be excluded a priori. These arguments can be ex- 
amined for the case of Hat-P-2b, for which the observed 
mass-radius relation requires a total mass of heavy ele- 
ments of at least 200 to 300 Mq, i.e. a mean mass- fraction 
Z> 7-11%. Hat-P-2b's parent star is a F8 metal-enriched 
{[Fe/H] = 0.1) star, with = 1.3 Mq. The maximum 
total amount of heavy element material available in the 
parent disk was thus about 900 Mq. So the required con- 
tent of heavy elements would be close to the aforementioned 
~30% upper limit of accretion efficiency. If both the mass 
and the density of Hat-P-2b are confirmed, this object thus 
lies at the edge of what is predicted to be possible within 
the current standard core-accretion scenario. We suggest, 
however, an alternative formation scenario, namely that 
the formation of Hat-P-2b involves collision(s) with one or 
several other massive planets. Besides forming big cores, 
collisions will lead to a substantial loss of the gaseous en- 
velope, thus to a larger relative fraction of heavy elements. 
Furthermore, gravitational scattering among planets gen- 
erally results in a tight orbit with a large eccentricity for 
one of the planets, which could explain Hat-P-2b's large 
eccentricity, and to the ejection of the other planet(s) or 
debris to interstellar medium. Such scattering processes be- 
tween planets seem to provide a viable and possibly dom- 
inant scenario to explain the observed eccentricity distri- 
bution of exoplanets ( [Chatterjee, Ford fc: Rasio 2007p . We 
thus speculate that Hat-P-2b was formed from such colli- 
sion processes. Note that a scenario based on giant impacts 
has also been suggested to explain the large heavy element 
content of HD149026b (Ikoma et al. 2006). 

Finally, in the present calculations, the heavy material is 
supposed to be composed entirely of one single component, 
water, rock (silicates) or iron. This is of course a simplifying 
assumption, as the inner composition of the planets is ex- 
pected to possess various fractions of each of these compo- 
nents. The water to silicate fraction, in particular, will vary 
depending whether the object has formed inside or beyond 
the ice line. Migration, however, will affect this fraction, 
as the migration process yields a larger collision rate of the 
planet embryo with rocky planetesimals and thus a decreas- 
ing abundance of volatiles. All ratios Mice/Mrock from to 
1 are thus probably possible. In any event, although the ex- 
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act composition of the heavy material component may have 
some impUcation on the mass-radius relationship for Earth- 
like planets (Valencia et al. 2006; Sotin et al. 2007; Seager 
et al. 2007), the present study shows that, for planets with a 
gaseous H/He envelope of mass fraction Mhhc > 10% Mp, 
current uncertainties on the EOS and the heavy element 
distribution lead to larger uncertainties on the planet's ra- 
dius determination than the effect due to different internal 
compositions (see Tables [T][2] and §5.1). 

7.2. The models 

In tables l4ll^. we present a subset of our grid of planetary 
models from 20 to 1 Mj with different levels of heavy 
element enrichment, Z = Zq, 10%, 50 % and 90%. For the 
purpose of the present paper, this grid is restricted for the 
moment to solar metallicity atmosphere models with two 
external atmospheric conditions, namely: (i) no irradiation 
(non-irradiated planets) and (ii) irradiation effects from a 
Sun at 0.045 AU, which is a typical incident irradiation for 
most of the transit planets discovered up to now. The effect 
of different levels of irradiation and different atmospheric 
compositions will be explored in a forthcoming paper. We 
have compared our models in this mass range with the mod- 
els of Fortney et al. (2007). We find an excellent agreement 
for jovian-mass planets using the same assumptions, i.e 
same core masses and level of irradiation. Small differences 
occur for Neptune-mass planets (Fortney's 17 model) 
for a large core mass fraction (Mcorc — 10 M®). This is 
due most likely to the zero-temperature EOS assumption 
adopted in Fortney et al. (2007) for the core EOS, as shown 
earlier in this paper (see §4.1). Larger discrepancies, how- 
ever, appear in the irradiated case for the Neptune-mass 
planet. This stems very likely from the different treatments 
of irradiation in the atmosphere models used in Forney et 
al. (2007, see Fortney & Marley 2007) and in the ones used 
presently (Barman et al. 2001), which seem to affect more 
drastically light planets. The exploration of these effects 
will be considered in a future paper. 

In the present grid of models, and in order to minimize 
the number of possibilities, the heavy element material is re- 
stricted to water using the SESAME EOS. The departures 
from this case due to different heavy material compositions 
(rock, iron), distributions within the planet's interior and 
EOS have been quantified in the previous sections of this 
paper (see Tables [T][2] for an illustration and §8 for a sum- 
mary) for the relevant planet-mass range. The model users 
are thus referred to these analysis to determine the varia- 
tions expected fom these various assumptions in the planet 
modelling. Our choice for the heavy element distribution 
in the model grid depends on the planet mass, and has 
been guided by the most realistic expected distribution, as 
discussed in §7.1 and summarized below, from our present 
analysis of the uncertainties resulting from different treat- 
ments of heavy element enrichment. The solar metallicity 
models, with Z — 0.02, have the same composition as that 
of a brown dwarf, but such objects with masses as low as a 
few Neptune-masses are of course not realistic. These mod- 
els are an extension to planetary masses of the models of 
Baraffe et al. (2003). Note that the 20 planet model 
with an enrichment as low a Z = 0.02 and with, conse- 

^ The complete grid of models, from 10 to 10 Mj, is avail- 
able on [http: / /perso. ens- lyon.fr /isabelle. baraffe /PL ANET08 , 



quently, a low mean density, expands under the effect of 
irradiation. This model is not included in Tab. [HI being 
meaningless. 

— For all masses and Z = 10%, all heavy elements are 
located in the core, since for such low Z, we have shown 
that their distribution has only a modest effect on the 
radius. 

— For larger enrichments {Z — 50%-90%), we make a dis- 
tinction between planets below Mp < 20 , hereafter 
denominated as " light planets" , and more massive plan- 
ets. Because all planets are expected to have a core of 
about ~ 10 Mg, the following distributions seem to be 
the most realistic: 

— For light planets, all heavy material is located in the 
core. 

— For more massive planets, we have shown that dis- 
tributing the heavy elements over the entire planet 
is similar to distributing them partly in a core of at 
most ~ 10% of the mass of the planet and partly 
in the envelope. We thus adopt such a distribution 
for these objects: the heavy material is distributed 
over the entire planet, using the AVL with SESAME 
EOS. The uncertainties due to the EOS (ANEOS 
versus SESAME) for this type of distribution have 
been quantified in §4 and §5. 

— For planets with masses > 1 Mj, we only provide models 
with Z=10%, since enrichments as large as 50% or 90% 
correspond to an amount of heavy material greater than 
the quantity available in protoplanetary disks around 
solar type stars. 

Figure shows the mass-radius relationships for plan- 
ets in the mass range 20 M^ - 1 Mj, for different levels 
of irradiation and heavy element enrichments. The models 
are compared to the observed mass-radius data of transiting 
planets. Note that in the case of irradiation (right panel), 
the theoretical planetary radius is not corrected from the 
effect of the atmospheric extension. Such a correction adds 
^ 4% to the calculated radius for this level of irradiation 
(Baraffe et al. 2003). 

8. Discussion and perspectives 

In this paper, we have explored the uncertainties in current 
planetary structure and evolutionary models arising from 
the treatment of the heavy material component in the in- 
terior of these planets. The study covers a mass range from 
10 M^ to 1 Mj. Our main results can be summarized as 
follows: 

— The ideal mixing entropy contribution arising from the 
heavy Z-material is found to be inconsequential on the 
planet's evolution, within the limit that the variation of 
the degree of ionization of these elements along the evo- 
lution is presently ignored. This mixing entropy, how- 
ever, represents about 10 to 20% of the total, H/Hc + Z 
elements, entropy, and thus modifies the internal isen- 
trope. A proper calculation of the planet's .structure for 
a given entropy should thus include this contribution. 

— For a metal mass fraction in the envelope Zcnv< 20%, 
the EOS of the Z-material can be approximated by using 
a corresponding Yoquiv effective helium mass fraction in 
the SCVH EOS. Above this limit, this approximation 
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Table 4. Radii of planets (in in the mass range 20 Mq - 
1 Mj for different levels of heavy element enrichment Z (see 
text, §7.2, for the choice of the heavy material distribution) 
and at different ages. 



100 200 300 
Mass [Mj] 



100 200 300 
Mass [Mj] 



Fig. 12. Mass-radius relationships for planets in the 
mass range 10 Mq - 1 Afj at different ages, as indi- 
cated in the panels, and different levels of heavy ele- 
ment enrichment: Z = Zq (solid lines); Z =10% (short- 
dash lines); Z =50% (dash-dot Hues); Z =90% (long- 
dash lines). Left panels: without irradiation; right pan- 
els: with irradiation from a Sun at 0.045 AU. Symbols 
indicate the observed data for transiting Neptune-mass 
to Jupiter mass planets, taken from F. Font's website 
(obswww.unige.ch/~pont/TRANSITS.htm). 



becomes more and more incorrect and yields erroneous 
cooling sequences. 

— For core mass fractions less than 50% of the planet's 
mass, a variation of the core composition between pure 
water and pure rock (iron) yields a difference on the 
radius of less than 7% (15%) after 1 Gyr, for all the 
planet masses of interest. 

— For a total mass fraction of heavy elements Z < 10%- 
15%, their impact on the evolution of the planet can be 
mimicked reasonably well by assuming that they are all 
located in the core. 

— For heavy material enrichments Z > 20%, the distri- 
bution of heavy elements (everything in the core versus 
uniform distribution) can affect significantly the cooling 
and thus the radius determination (more than 10% at 
a given age). Therefore: 

— For metal-rich {Z > 20%) light planets {< 20 M(j)), since 
the planets are expected to have a massive ~ 10 M0 
core, it seems realistic to put all the heavy material in 
the core. 

— For massive metal-rich planets (Mp > 50 and Z > 
20%), however, the evolution is better described by 



z 


Mp/Me 


-R0.5 Gyr 


Rl Gyr 


R5 Gyr 


0.02 


10. 


0.828 


0.811 


0.758 




20. 


0.858 


0.839 


0.800 




50. 


0.923 


0.905 


0.876 




100. 


0.980 


0.963 


0.937 




159. 


1.017 


0.995 


0.968 




318. 


1.057 


1.032 


0.998 


0.10 


10. 


0.779 


0.763 


0.716 




20. 


0.813 


0.797 


0.762 




50. 


0.878 


0.862 


0.836 




100. 


0.935 


0.919 


0.896 




159. 


0.971 


0.951 


0.926 




318. 


1.012 


0.990 


0.958 


0.50 


10. 


0.598 


0.586 


0.555 




20. 


0.632 


0.621 


0.598 




50. 


0.683 


0.665 


0.635 




100. 


0.717 


0.697 


0.670 




159. 


0.739 


0.718 


0.690 




318. 


0.781 


0.756 


0.721 


0.90 


10. 


0.382 


0.375 


0.357 




20. 


0.420 


0.414 


0.403 




50. 


0.503 


0.474 


0.439 




100. 


0.543 


0.512 


0.469 




159. 


0.568 


0.538 


0.492 




318. 


0.607 


0.578 


0.524 


Table 


5. Same 


as table |4] for irradiated mc 


at 0.045 AU. 








Z 


Mp/Me 


^^0.5 Gyr 


Rl Gyr 


R5 Gyr 


0.02 


10. 










20. 


1.441 


1.391 


1.229 




50. 


1.258 


1.192 


1.084 




100. 


1.201 


1.155 


1.074 




159. 


1.186 


1.151 


1.085 




318. 


1.160 


1.137 


1.089 


0.10 


10. 


1.517 


1.506 


1.428 




20. 


1.416 


1.326 


1.147 




50. 


1.173 


1.120 


1.021 




100. 


1.131 


1.091 


1.020 




159. 


1.120 


1.089 


1.032 




318. 


1.106 


1.086 


1.042 


0.50 


10. 


1.095 


1.050 


0.927 




20. 


0.935 


0.890 


0.800 




50. 


0.841 


0.798 


0.725 




100. 


0.812 


0.787 


0.737 




159. 


0.813 


0.790 


0.747 




318. 


0.830 


0.810 


0.772 


0.90 


10. 


0.545 


0.528 


0.492 




20. 


0.517 


0.502 


0.474 




50. 


0.575 


0.534 


0.469 




100. 


0.594 


0.560 


0.498 




159. 


0.601 


0.574 


0.522 




318. 


0.628 


0.604 


0.552 



models which assume that the heavy elements are dis- 
tributed throughout the entire planet than by models 
with all heavy elements in the core and none in the 
gaseous envelope. The former models yield results sim- 
ilar to the ones obtained with a more realistic distri- 
bution, namely a ~ 10 core and the rest of heavy 
material distributed in the envelope. 
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— The tcmpc;r;iture dependence of the heavy material EOS 
and the release of gravitational and thermal energy of 
the core have negligible effects on the cooling history 
of m,assive planets (saturnian and jovian masses), inde- 
pendently of the core mass. For Neptune mass planets, 
these effects are significant (~ 10% difference on the ra- 
dius after 1 Gyr) for extreme heavy element enrichments 
[Z > 90%), as for instance in GJ 436b. 

We have performed comparisons between available EOS 
of heavy material which cover a range of temperature and 
pressure large enough to cover the characteristic domains 
of exoplanet interiors. Unfortunately, significant differences 
exist between these EOS concerning the thermodynami- 
cal quantities relevant for planet evolution, such as the en- 
tropy, in particular at high pressures (P 0.01 Mbar) and 
temperatures (T ^ 5000K), where no experimental or nu- 
merical guidance is provided. These discrepancies translate 
into non-negligible variations on the mass-radius relation. 
Therefore, by venturing into parts of the phase diagram of 
heavy elements that are presently not possible to reach ex- 
perimentally or with computer numerical simulations, the 
planet internal structure calculations include inevitably a 
certain degree of uncertainty, which has been quantified as 
thoroughly as possible in the present study. For planets 
with even a modest fraction of H/He gaseous envelope, as 
expected above '-^6 — 10 M©, these heavy material EOS 
related uncertainties are the major culprit for preventing 
accurate determinations of the exoplanet internal composi- 
tion from the observed mass and radius. "Accurate" means 
in this context at a level better than yielding < 10% varia- 
tions on the radius. For Earth or super-Earth planets with 
no gaseous envelope, a radius measurement accuracy bet- 
ter than 5% is expected to allow to distinguish icy from 
rocky internal compositions (Valencia et al. 2007). However, 
as shown in §5 for the specific cases of HD 149026b and 
GJ436b, in the presence of even a small (< 30% by mass) 
H/He contribution, varying the distribution of the heavy 
material within the planet has by itself a larger impact, 
not mentioning the one due to the uncertainty in the EOS. 
Therefore, for planets above ~ 10 M©, massive enough to 
accrete an H/He envelope, it seems difficult to determine 
precisely the internal composition with current structure 
models, as the effect of the heavy material composition on 
the radius is blurred by the presence of a gaseous enve- 
lope. Even a modest fraction (~ 10%) of H/He already 
severely modifies the evolution and thus the radius of a 
planet compared with a gasless one. As an example, a 10 
M© mass planet retaining a 10% H/He envelope is ~ 1.5 
times larger than its pure icy counterpart (Seagcr et al. 
2007). Frustratingly enough, the impact of this gas contri- 
bution casts doubt on our ability to determine accurately 
these planet inner composition from the observed radius, at 
a better level than figuring out their gross bulk properties, 
already an important information for our understanding of 
planet formation. 

The main goal of this work was to explore thoroughly 
current uncertainties and assumptions in current models 
of planet structure and evolution. Such a study had been 
eluded so far in the field of exoplanets and becomes now 
necessary, given the high accuracy achievable by present 
and forthcoming observations. This work enters and pro- 
longates the long history of the modeling of giant planets, 
where some of the problems and uncertainties discussed 



in the present paper have been addressed decades ago for 
the case of our own Solar System planets. As an example, 
the sensitivity of Uranus and Neptune models to the EOS 
of water at high pressure has been thoroughly discussed 
by Hubbard & MacFarlane (1980). The importance of the 
thermal and gravitationnal energy release of heavy material 
in planet interiors was already stressed for the evolution of 
Uranus and Netpune (Podolak et al. 1991; Hubbard et al. 
1995). All this historical work has already highlighted the 
remaining large uncertainties on the modeling of our own 
giant planets, despite a wealth of observational constraints 
which are lacking for exoplanets. 

Although our conclusions sound rather pessimistic, con- 
cerning the degree of accuracy expected from current planet 
models, this work should motivate further efforts both on 
the experimental and theoretical fronts, in order to make 
progress in this thriving field of exoplanet exploration. The 
(non exhaustive) homework list of improvements includes 
the following items, (i) Exploring the EOS of heavy ma- 
terial in the critical pressure regime 0.1-100 Mbar and at 
high temperature is crucial. Substantial progress in this do- 
main is expected with the ongoing and future high-pressure 
experiments in various national laboratories (including e.g. 
the Lawrence Livermore and Sandia labs in the US or the 
LIL and Mega Joule laser projects in France), (ii) In the 
same vein, first-principle N-body numerical methods (DFT, 
path-integral, quantum molecular dynamics, ...) should be 
able to provide at least some benchmarks in the part of the 
P-T diagram which for now lies in the unknown interpo- 
lated regime, (iii) It is important to explore the efficiency of 
the heat transport mechanism in planetary interiors in the 
presence of molecular weight gradients, in particular in the 
case of large heavy material enrichment. This bears crucial 
consequences on the planet's cooling history (Chabrier & 
Baraffe 2007). Current progress in high-resolution multi-D 
numerical simulations should be able to handle this prob- 
lem, (iv) Even though measuring the extrasolar planet grav- 
itational moments is hardly conceivable, it might be pos- 
sible to determine their oblateness, e = Req/{Req — Rpoi), 
where Req and Rpoi denote the equatorial and polar ra- 
dius, respectively, with future transit observations (Seager 
& Hui 2002; Barnes & Fortney 2003). This in turn leads to 
the determination of the rotation rate and thus of the cen- 
trifugal potential, providing a more stringent constraint on 
the internal structure of the planet, (v) Last but not least, 
the perspective, on the observational front, of direct plan- 
etary atmosphere observations (LYOT project, GEMINI, 
ELT, DARWIN/TPF or their precursors) and transit detec- 
tions (CoRoT, Kepler) will improve our knowledge of their 
surface composition and radius measurements and provide 
important constraints on the planet's content in heavy ma- 
terial. 

Finally, we have shown in this paper that massive ( ^ 
several Jupiter mass) planets may form from two differ- 
ent avenues, namely the standard core accretion scenario 
(Mordasini et al 2008), and giant impacts between massive 
planets or planet embryos. This latter process is likely to 
yield very metal-enriched and thus very dense massive plan- 
ets, with a finite eccentricity, as a result of planet scatter- 
ing. Objects like HD 149026b; Hat-P-2b or the very recently 
discovered HD 17156b (Gillon et al. 2007c) could be the il- 
lustrations of this latter planet formation mechanism. The 
observation of Hat-P-2b, together with the mmierous obser- 
vations of free fioating brown dwarfs of a few Jupiter-masses 
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(Caballero et al. 2007) shows that planets and brown dwarfs 
have a substantial (about one order of magnitude in mass) 
mass domain overlap. As we have shown in this paper, 
planets massive enough to exceed the deuterium-burning 
mass limit will indeed ignite this reaction at the bottom 
of their H/He rich envelope, at the top of the core. This 
is one more evidence, if it were still necessary, that using 
the deuterium-burning limit as a criterion to distinguish 
planets from brown dwarfs has no valid foundation. 

The complete grid of models, from 
10 Mffl to 10 Mj, is a vailable on 
http://perso.ens-lyon.fr/isabelle.baraffe/PLANET08i 
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Appendix A: Calculation of the mixing entropy 

In this appendix, we calculate the expression for the ideal 
mixing entropy of a two-component system composed of a 
H/He mixture, identified as one component, on one side and 
of a heavy material component of mass fraction Z — Mz /M 
on the other side. In the following, the subscript " 1" denotes 
the H/He component while the subscript "2" refers to the 
Z-component. By definition, the ideal entropy of mixing 
reads: 



^ M\nN - NiXiyNx - N2\nN2 - A^elniVe, (A.l) 

Kb 

where M = Ni + N2 denotes the total number of parti- 
cles, including free electrons, in the fluid, J\fi = Ni -\- Nl 
denotes the total number of particles of component z, with 
Ni the number of nuclei of component i and Nl the number 
of electrons provided by the component i. Developing eq. 
(|A.ip yields: 

=Mln(l-f^)-fAA2ln(l-f^) 
+ MlnWi +A/'2lnA/'2 

- iVi IniVi - N2 IniVa - N^ InTVe 
= Mln(l + ^)+AA2ln(l + ^) 

- Ne\iiNe + NllnNl 

+ Nl\nNl+'-f^+'-f^ 



20 



Baraffe, Chabrier & Barman: From super Earth to super Jupiter exoplanets. 



where 



: M In M - N, In N, - Nl In Nl 



(A.2) 



denotes the ideal mixing entropy of the component z, in- 
cluding various ionic, atomic or molecular species as well 
as electron contributions (see SCVH) , and Ne ~ Nl + N"^ 
is the total number of electrons provided by the H/He and Z 
components. These ideal mixing entropy contributions are 
already included in the SCVH EOS for the H/He compo- 
nent (their eq.(53)) Q and in the appropriate EOS for the 
Z-component. Removing these two contributions, we ob- 
tain the ideal entropy of mixing which arises only from the 
mixture of the H/He and Z components: 



Nl ln(l 



f ^)-fA/-2ln(H-^) 
+ Nl\^N} + N'i\^N^ 



(A.3) 



The specific entropy, i.e. the entropy per unit mass is 
given by 

kB " MkB " VlUkB ^ 

where M — Mi + M2 is the total mass. 

We define a mean atomic mass = Mi/Ni and mean 
charge 1 + Zi = Ni/Ni for each component, where we have 
used the electroneutrality condition Nl = NiZi, so that 
^ = rni/{l + Zi). For the H/He component, the quanti- 
ties mi and Zi are given in terms of the relative number 
fractions a; of H2, H, He, He+, He++ and e " by the SCVH 
EOS (then eqs.(33)-(35)), with: 



After some algebra, the ideal specific entropy of : 
of the (HHe)/Z mixture can be written: 



MkE 



|^(l + Zi)x <j ln(l-H/3Y)+/3Yln(l + ^) 



xi'hn(l + <5')-/3'7'^e'^Ml + i7) 



(2) 



where Xi = Mi/M = Mhhc/M and 



/3' 



mi Z 
m2 1 - J 



I + Z2 
l + Zi 



with Z = Mz/M and x'i^ = Z,/{l + Z., 



m 



I Z2 J^e „/ , 

(A.8) 



„(2) 



Ml , , , 

—J- = (2a; ^2 -I- a;^ -I- xh+ jmn + [XHe^+ + XHe+ + XHejrriHe 
All 

(A.5) 
(A.6) 



and 



„(i) 



Zi = 



1 — a;, 



(1) 



where 



X 



(1) 



Nl 1 , /37 



Ui I + /37 " l + /?7 



XH++ , , o ixHe++2xHe2+) (A. 7) 



denotes the number-concentration of free electrons in the 
HHe mixture as defined in SCVH, niH and rriHe denote the 
atomic mass of hydrogen and helium, respectively, and the 
grec symbols a,/?, 7 are defined by eqs.(54)-(56) of SCVH. 
For the Z-component, the mean mass corresponds to the 
atomic mass of the compound under consideration. The 
mean charge Z2 is unknown. We have carried out calcula- 
tions for the two limiting cases of fully neutral {Z2 = 0) 
and fully ionised (Z2 = Znuc) heavy material, where Znuc 
denotes the nuclear charge of the compound. 



Note the following typos in the SCVH paper. In eqs. (45)- 
(46), the entropy ratios on the rhs of the equations should be 

^ and ^-g- and not and -gHF- lu eq.(56) for the parameter 
5, the fraction in front of the bracketed terms on the rhs of the 
equation should be | and not |. 



